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Abstract 

In this comment, I point out to several mathematical mistakes in the above-referenced 
letter. 



In a recent letter (Jl)), Xu, Huang and Yu (referred to as XHY below) have 
derived effective nonlinear susceptibilities of graded composite films. The the- 
oretical approach used by XHY is based on the perturbative theory developed 
previously in Ref. (0) by Yu (the same author as in the letter which is subject 
of this comment), Hui and Stroud. However, XHY make several mathemati- 
cal mistakes in applying the formalism of Ref. pi), which are briefly detailed 
below. 

XHY have considered a two-component film whose composition is varied in 
one dimension. More specifically, the volume fractions of the two components 
making up the film, p\ and p 2 (pi + P2 = 1) were assumed to depend on the 
linear coordinate z G [0, L\. NeverthelesSjthe computation of effective nonlin- 
ear responses, given in Eqs. 1-15 of Ref. (jlj), was carried out for fixed values of 
Pi and p2- The dependence of effective linear and non-linear susceptibilities on 
z was then expressed through the dependence of p\ and p 2 on z, as is evident 
from Eqs. 16,17 of Ref. (0). Thus, the gradation of the films was assumed to 
be slow enough, so that the effective susceptibilities could be physically de- 
fined as functions of z. This assumption can be reasonable and is not subject 
of this comment. However, the method used for the derivation of the effec- 
tive constants at fixed values of pi and p 2 outlined in Eqs. 1-15 of Ref. (0) is 
erroneous. 

The derivation mentioned above is based on Eq. 14 of Ref. (0). This equa- 
tion gives a formula for computing the effective linear dielectric constant of a 
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composite, e e , and, written in terms of the volume fractions p\ and p 2 and the 
linear dielectric constants €\ and e 2 of the components, has the form 

{El) (E 2 2 ) 



In this formula, (Ef) and (E%) are the averages of the square of the (real- 
valued) linear electric field, computed inside the first and the second con- 
stituent of the film, respectively, and Eq is the external electric field. Because 
Ei and E 2 are linear fields, the ratios (Ei 2 )/Eq are independent of E . Thus, 
Eq. 1 gives the effective linear dielectric constant of the composite in terms of 
the linear dielectric constants of its constituents as a weighted average. Inci- 
dentally, this definition of e e was obtained by equating the total electrostatic 
energy of a homogeneous sample occupying some volume V and character- 
ized by the linear dielectric constant e e and that of a composite occupying 
the same volume and characterized by the constants p 1)2 and ei j2 . Of course, 
the ratios (E^ 2 ) / Eq must be computed by solving the electrostatic boundary 
value problem for each specific geometry of the composite. 

In what follows, I detail three mistakes XHY have made in applying Eq. 1 to 
the the problem of computing the nonlinear responses of composite films. 



The First Mistake At the onset, XHY use Eq. 1 to compute the av- 
erages (El) and (E 2 ), assuming that e e is given by some known function, 
e e = F(ei, e 2 ,pi,P2)- Namely, they write (Eqs. 2,3,7 of Ref. |l|)): 

(El) = l dF &>*»Pi>P*) E * ( 2 ) 

Pi OEi 



and analogously for (E 2 ). In writing this equation, XHY have ignored the fact 
that the averages (Ef) and (E 2 ) are themselves functions of both ei and e 2 . 
Differentiation of Eq. 1 with respect to e\ (taking into account e e = F) would 
yield 

n (lF*\ I r d{El) \ I „ r _ dF(e 1 , e 2 , Pl , p 2 ) 

Pl [(Ex) + ei-^ r ) +^ 2 — ^ E . (3) 

Eq. 2 is derivable from Eq. 3 only if d(Ef)/dei = d(E 2 )/dei = 0, which, 
obviously, is not the case. 



The Second Mistake Application of Eq. 2 requires the knowledge of the 
function F(ei, e 2 ,pi,p 2 ). To this end, XHY define F for the graded film as 
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a whole (not locally) by writing F = L[Jq e^ G (z)dz]~ 1 (Eq. 21 of Ref. (0)), 
where euc is the Maxwell-Garnett effective dielectric constant that can be 
found analytically from Eqs. 16,17 of Ref. (Q). The first step in this procedure 
has not been justified in Ref. (jl|) and appears to be arbitrary. More impor- 
tantly, the second step requires that the Maxwell-Garnett formula give the 
same result for the dielectric constant as Eq. 1 with properly computed field 
averages. But the two equations are, generally, not equivalent. In Ref. (@), a 
Maxwell-Garnett-type formula was derived from Eq. 1 for the case of spheri- 
cal inclusions of volume fraction pi in the limit p\ — > 0. But for mixing ratios 
close to 0.5, and for high-quality metal inclusions in a dielectric host with the 
electromagnetic frequency being close to the Frohlich resonance of a single 
inclusion (all of which is the case in numerical examples shown in Ref. (JlJ)), 
the Maxwell-Garnett formula is known to be very inaccurate 0). In fac"it 
has been demonstrated in Ref. (0) that the Maxwell-Garnett theory provides 
a reasonable approximation only for p\ < 10 -3 . 

Physically, computing the averages (Ef) and {E 2 ) by means of Eq. 2 (which is, 
also incorrect due to the First Mistake), where the analytical form of F is de- 
rived from the Maxwell-Garnett formula ignores the well-known phenomenon 
of strong fluctuations of electric field in resonant composites (0; M) ■ 



The Third Mistake XHY make the most serious mistake when they state 
that the nonlinear susceptibilities can be obtained from Eq. 1 by viewing 
the constants as intensity-dependent, i.e., by making the substitutions 
ei,2 - > ei,2 + Xi,2(Ef 2 ) i n t ne arguments of the function F, and by expanding 
F with respect to the small parameters Xi,2- This procedure is mathematically 
expressed in Eq. 6 of Ref. (jlj) and is clearly erroneous. By definition, ei j2 in 
Eq. 1 are the linear dielectric constants of the composite constituents, as well 
as e e is, by definition, the linear effective dielectric constant of the composite. 
Computation of higher-order effective susceptibilities would require computing 
higher moments of the field. For example, the formula derived in Ref. (0) for 
the effective third-order susceptibility is 

(Et) (£ 2 4 > 

Xe =PlXl L ^ L +P2X2^ ■ (4) 



Thus, Eq. 6 of Ref. (|l|) is based on an incorrect interpretation of the theoretical 
results of Ref. (0) 

Finally, the conclusion of Ref. (|lj) that "the harmonics [generated in a graded 
composite film] are significantly dependent on the gradation profiles as well as 
the width of the composite film" is technically (and trivially) correct. However, 
equations derived in Ref. (j[l) can not be used to "monitor the gradation profile 
as well as the width of the composite graded film by measuring the nonlinear 
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ac responses of the film subjected to an ac electric field" due to the errors 
described in this comment. 
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